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Abstract 


This  study  was  sponsored  by  DARPA,  in  the  wake  of  last  year’s  ONR-sponsored 
JASON  Study  JSR-99-215  on  Fast  Transoceanic  Transport  (aka,  “Fast  Ships”  study). 

The  charge  of  this  year’s  study  was  to  focus  on  turbulent  boundary-layer  drag  reduction, 
which  was  identified  as  necessary  for  transoceanic  transport  at  high  speeds,  i.e.,  at 
speeds  U>75  knots. 
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Previous  JASON  Work  —  / 


PED  et  al.  1999  “Fast  Ships:  Hydrodynamics  of  Fast  Ocean  Transport  ”  JSR- 
99-215. 


-  Feasibility  study  of  transport  vessels  characterized  by 

*  high-speed:  75  knot  + 

*  long-range:  10  knm 

*  Moderate  displacement:  10  ktons 


>  This  requires  (system)  LID  »  40 

*  LID  «  1 8-20  achievable  today 

*  Roughly  70%  of  drag  attributable 
to  boundary-layer  friction 

o 

*  LID  «  40  requires  «  x3-4  2 

reduction  in  BL  drag. 


JASON  2000  Drag  Reduction 


The  previous  JASON  Study  (JSR-99-215)  on  fast  transoceanic  transport  considered 
a  high-span  (6  —  65  m),  dual-hydrofoil-bome  ship  as  a  notional  design  as  perhaps 
the  only  viable  choice  at  the  speeds  indicated.  It  concluded  that  a  system  lift-to- 
drag  ratio  of  LID  «  40  was  necessary  to  approach  the  desired  specifications  of  a  fast 
Army  transport. 

Economic-viability  issues,  e.g.,  that  a  commercial  carrier,  such  as  FedEx,  should  be 
able  to  operate  a  fleet  of  such  vessels  also  argue  for  a  high  -LID  vessel.  That  was 
also  a  consideration  last  year’s  study  was  invited  to  consider. 

The  figure  plots  LID,  against  percentage  skin-friction  drag  reduction,  for  a  ship  that 
achieves  LID  =  20  without  drag  reduction,  assuming  that  the  unreduced  skin  friction 
is  70%  of  the  total  drag. 

LID  *  20,  achievable  with  high  optimization  of  conventional  technology,  is  still  an 
ambitious  goal.  The  LID  for  a  Boeing  747  —  a  very  capable  flyer  —  is  in  the  range 
of  16-18  and  included  as  a  reference. 
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What  Can  Drag  Reduction  Do  For  You?  —  Range 


The  range  relation  plotted  derives  from  the  Breguet  equation  ( e.g .,  JSR-99-215) 

R 

n.mi 

in  which, 

r)p  «  0.7  is  the  propulsion  efficiency  (plant  HP  to  propulsive  HP), 

LID  is  the  (average)  lift-to-drag  ratio  (treated  as  a  constant), 
wfucl  is  the  initial  displacement  (weight)  fraction  allocated  to  fuel, 
and 

cfucl  «  0.4  lbm/HP-hr  is  the  (reciprocal)  of  the  fuel  energy  content  per  unit  mass. 


v\p  LI D 

cfUCi  /  [lb  /  HP  •  hr] 


ln  (1_Wfuel) 
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The  expression  below  the  figure  and  the  corresponding  plots  express  the  required 
(plant)  power  (force  x  speed),  including  the  consequence  of  the  propulsion 
efficiency,  %,  in  kHP, 

If  high  speed  is  to  be  attained,  it  is  the  power  requirements  that  dictate  the  high  L!D 
values  mentioned  above  and  preclude,  as  a  practical  matter,  contemplating  speeds  ii 
excess  of  75  knots. 


Previous  JASON  Work  —  II 


•  P.  Diamond  et  at.  1992  “Drag  Reduction  by  Polymer  Additives,”  JSR- 
89-720. 

>•  Extending  Bird  et  at.  (1987),  proposed  a  coupled  tensor- field  description 
of  the  local  elastic-elongation  response  of  a  dilute  polymer  solution  to 

•  hydrodynamic  stretch  (strain  rate), 
opposed  by 

•  the  characteristic  relaxation  time  of  stretched  polymer  “dumb-bells” 

>  Noted  formal  analogy  between  coupled  hydrodynamic  and  polymer- 
stretch  fields  with  MHD  flows. 

-  Did  not  address  drag-reduction  mechanism(s). 

•  C.  Callan  and  K.  Case  1981  “Drag  Reduction  and  Adaptive  Boundary 
Conditions,”  JSR-8 1-17. 

-  Boundary-layer  stability  and  potential  increases  of  stability  margins. 


J A  SO  IS  2000  Drag  Reduction 


The  Diamond  et  at.  (1992)  JSR-89-720  report  provided  an  important  summary  of 
ideas  at  the  time,  as  well  as  a  significant  impetus  and  guidance  for  this  year’s  study. 
It  identifies  the  mechanisms  and  dynamics  that  couple  polymer-elastization  fields 
and  turbulence. 


The  Callan  &  Case  (1981)  JSR-81-17  study  focused  on  hydrodynamic  stability  and 
transition  out  of  the  laminar  state  and,  as  a  consequence,  is  less  germane  to  the  high¬ 
speed  regime  of  interest  in  the  present  context. 
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Three-  to  Four-Fold  Drag-Reduction? 


•  Technologies  with  demonstrated  x3-4  drag-reduction  potential 

-  High-MW  polymers 

-  Surfactants  (surface-active  agents) 

-  Microbubbles  (dh  <  50  pm) 

-  Air/vapor  films 

•  Other  technologies  considered 

-  Electromagnetic  Turbulence  Control  (EMTC) 

-  Active  (feedback)  control 

•  Issues 

-  Momentum  transport 

-  Drag-reduction-power  expended  vs.  propulsive-power  saved 


JASON 2000  Drag  Reduction 


In  addition  to  the  four  technologies  listed  at  the  top  that  meet  the  x3-4  drag- 
reduction  criterion,  two  other  technologies  were  considered. 

Electromagnetic  Turbulence  Control  (EMTC)  was  cited  by  Du  &  Kamiadakis 
(2000,  Science)  in  their  computer  simulations  of  moderate-i?e  boundary-layer  flows 
as  capable  reducing  drag  by  factors  comparable  to  the  values  of  interest  here,  with 
small  required  powers.  John  Kim,  in  his  July  2000  briefing  to  us  (also  Kim  2000), 
agreed  with  the  drag-reduction  findings  by  Du  &  Kamiadakis,  but  offered  a 
differing  estimate  of  power  requirements  for  a  given  drag-reduction  level. 

Resolving  this  disagreement  was  one  of  the  issues  examined  as  part  of  this  year’s 
study. 

Active  (feedback)  control  technology  has  also  recently  claimed  large  drag-reduction 
benefits,  based  on  computer  simulations. 

Issues  addressed  included  mechanisms  for  momentum  transport  and  turbulent 
boundary-layer  drag  reduction,  as  well  as  issues  of  efficiency.  Efficiency  here 
means  the  ratio  of  power  saved  by  drag  reduction  to  power  expended  to  achieve  it. 
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Momentum  For  A  Continuum  Material 


Momentum  conservation  for  a  differential  volume  element  of  a 
continuum  substance  (Cauchy’s  equation) 

Du 

p -  =  pf-Vp+V-T 

y  Dt 

where, 

-  p  =  p(x,f)  is  the  density  field, 

-  u  =  u(x,r)  is  the  velocity  field, 

-  D/DfH=d/df  +  u-  V  (convective  derivative), 

-  f  =  f  (x,/)  is  any  body-force  field  ( e.g .,  E&M  Lorentz  force), 

-  p  =p(x,t)  the  pressure  field,  and 

-  t  =  t(x,/)  the  total  stress  tensor  field 

•  x ij  ~  Zji  for  isotropic  fluid 


|  JASON  2000  Drag  Reduction _ | 

In  considering  drag-  and  momentum-transport  issues,  the  general  (Cauchy)  equation 
should  be  used,  in  which  the  form  of  the  body  forces  and/or  stress  tensor  can  be 
defined  as  appropriate  to  capture  imposed  external  body  forces  as  well  as 
modifications  to  the  flow  stress  field  as  a  consequence  changes  in  the  rheology  of 
the  fluid.  In  this  form,  it  applies  to  the  motion  of  both  Newtonian  and  non- 
Newtonian  fluids. 

D/D?  denotes  the  convective  (Lagrangian)  derivative  and  be  expressed  in  terms  of 
local  (Eulerian)  time  and  space  derivatives,  i.e., 


■  +  u-V 


As  expressed  in  the  equation  in  the  frame,  f  is  a  body  force  field  and  has  units  of 
acceleration,  e.g.,  gravity. 

For  an  isotropic  fluid,  the  stress  tensor  must  be  symmetric.  Isotropy  may  or  may  not 
be  valid  as  an  assumption  for  a  solution  of  long-stranded  molecules,  such  as  high 
molecular- weight  polymers. 
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Stress  Tensor  —  Newtonian  Fluid 


For  a  Newtonian  fluid,  the  only  stresses  are  viscous,  i.e.. 


t  =  -r(visc)  - 


=  2F5 


where, 


D 


+  H,(V-u)l 


jis  is  the  (shear)  dynamic  viscosity, 

-  |iv  is  the  volume  (bulk)  viscosity  (viz.  sound  attenuation), 

-  D  =!4[Vu  +  (Vu)T  ],  is  the  (rate  of)  deformation  tensor, 

-  I  is  the  identity  tensor,  and 

-  V  u  ==  tr{D}  is  the  velocity  divergence,  where  (mass  conservation), 


1  Dp 

- -  =  -V  u 

p  Dt 


JASON  2000  Drag  Reduction 


This  expression  for  the  stress,  a  function  of  two  viscosity  coefficients,  defines  a 
simple  (single-component)  Newtonian  fluid.  The  shear  and  bulk  viscosity 
coefficients,  ps  and  pv  [Pa-s  in  MKS],  are  thermodynamic  variables,  and  functions 
of  the  fluid  and  not  the  flow. 

In  Cartesian-component  form,  the  deformation  tensor,  D,  and  velocity  divergence 
(dilatation),  V-u,  are  given  by, 


KdxJ 


+*C| 


dxi  J 


V  n 


In  the  context  of  Fast  Ships,  including  the  contribution  of  dilatation  (V-u)  is 
necessary  in  the  study  of  microbubbles,  which  increase  compressibility  and 
decrease  the  speed  of  sound  substantially,  leading  to  significant  compressibility 
effects  at  the  speeds  of  interest  ( U  >  35  m/s).  This  will  be  discussed  later. 
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Incompressible  Newtonian  Fluid 

For  constant-density  flow  (V  u  =  0)  and 
T(visc)  =  2|isD  . 

for  f  =  0  and  ps  *  fn(x),  with  v  =  ps  /  p,  this  yields  the  familiar, 
innocent-looking  equation, 

—  =  ~-yp  +  vV2u 
D  t  P 

which,  as  far  as  we  know,  describes  (all  of)  incompressible  turbulence 
for  a  Newtonian  fluid. 


JASON  2000  Drag  Reduction 


In  the  context  of  Fast  Ships,  for  drag-reduction  not  based  on  microbubbles  or 
air/vapor  films  endowed  with  significant  dynamic-pressure  variations,  dilatation  is 
negligible  and  the  incompressible-flow  equations  can  be  used. 
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For  laminar  boundary  layers  (Blasius-equations  solutions),  with  x  the  wetted  length 
(upstream  streamwise  extent  of  growing  boundary  layer),  we  have: 


tW‘- 


1 721  6 

1/2  =  0.34—  (displacement  thickness) 


^ ,  iUi-i  . 

V  V  j  '  u{  U  I 


0.664  A  .  8  , 

— j -pj-  =  0.13—  (momentum  thickness) 


pf/2/2 


8,  0.664  , , .  . .  .  ^  N 

=  —  =  (skm  -  friction  coefficient) 


Cd(L)  S  Jdx  Cf(x) 


Re,  = 


0 

p  U  x 


1.328 

Re l'2 
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Boundary-Layers  —  Turbulent  Flow 


treestream  velocity  - 

u(y) 

y* 

1  u 

s  velocity 

bound  aiy-1  ay  er  edge - n. 

- J 

profile 

5 

wall 


Separate  into  mean  and 
time-fluctuating  quantities 

p(x,/)->p(x)  +  p'(x,/), 
u(x,0->u(x)  +  u'(x,0, 

A 


Boundary-layer  approximation  (Reynolds-averaged  equations): 


du  +  dv  _  ^ 
dx  dy 

du  du  TTdU  1  dxcfr 
dx  dy  dx  p  dy 

p(x)  =  p(x)-py'2 


For  a  Newtonian  fluid, 


du 


with, 


*efr=  P—  ~  pwv  . 
dy 

UL  ’ 


Reynolds 

stress 


^wall  P 


w(x,0)  =  v(x,0)  =  0 ,  u(x ,  oo)  =  JJ .  same  as  for  a  laminar  boundary  layer, 

since  u'  -»  0  and  v'  0,  as  y  ->  0. 

JASON  2000  Drag  Reduction 


Phenomenologically,  in  terms  of  their  action  on  mean  flow  profiles,  the  cross¬ 
stream  correlations  of  the  turbulent  velocity  fluctuations  can  be  viewed  as 
contributing  an  additional  stress,  dubbed  the  Reynolds  stress. 

For  flow  of  a  Newtonian  fluid,  the  shear-stress  at  the  wall,  tw,  has  the  same  form  for 
steady  (laminar)  and  turbulent  flow,  because  fluctuations  vanish  as  the  wall  is 
approached,  driving  Reynolds  stresses  to  zero. 
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Boundary-Layers  —  Velocity  Profiles 


*  Coupling  of  xw  to  velocity  profile  applies  to  Newtonian  and  non-Newtonian 
fluid  flow,  if  other  stresses  vanish  at  the  wall 

*  Changes  in  shear-stress  at  the  wall  linked  to  changes  in  boundary-layer 
velocity  profile. 

—  Potentially,  a  rock  and  a  hard  place:  in  decreasing  tw  one  can  produce  unstable  BL 
velocity  profiles 

*  An  integral  relation  between  BL  velocity  profiles  and  shear  stress  at  the  wall 
can  be  expressed,  for  general  fluids,  by  extending  the  T.  von  Karman 
boundary-layer  momentum  integral. 

JASON 2000  Drag  Reduction 


For  Resi  >  104,  or  so,  II »  0.62.  For  a  smooth  wall,  the  skin-friction  coefficient  can 
be  approximated  by  the  Prandtl-Karman  relation. 


a  4.0  log 5.6 


or,  the  explicit,  modified  (JSR-99-215)  Hama  (1954)  relation, 
1 


4$ 


=  3.44  log10i?ex  -  5.6 


Turbulent  boundary-layer  thickness  scales  are  given  by, 


8 

x 


h. 

8 

h. 

8 

A 

82 


0.16 

Reln 


a 


\Uj 
\ 


=  a^/Cf /2 


a 


ft 

M 


U 


l-6.8^Cf/2 


a  =  3.6; 
p  =  6.8  a; 
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Boundary-Layer  Momentum  Integral 


Can  integrate  BL  equations  to  generalize  the  von  Karman  (1921)  result 
(p  here  assumed  uniform) 


pc/2 


2 


dS2  _1_ 
dx  +U 


(282  +  5,)^-vw 


1 

U 2 


\fAy 


-  For  both  Newtonian  and  non-Newtonian  fluid  flow 

-  U-  U(x )  is  the  freestream  velocity  (df/ /dx  =  0  if  dp/dx  =  0), 


dy 


and 


dy 


-  Cf  is  the  skin-friction  coefficient, 

-  8j  and  S2  are  the  displacement  and  momentum  thicknesses  [  =  0(8)], 

-  vw  is  the  wall-normal  local  blowing/suction  velocity  (+  for  injection), 

-  fx  is  the  body-force  contribution  (e.g.,  gravity  or  Lorentz  force). 

•  Must  be  included  in  estimating  tw  in  the  presence  of  MHD  forces 

•  Indistinguishable  from  MHD  propulsion. 


JASON 2000  Drag  Reduction 


The  result  derived  is,  basically,  the  same  as  the  von  Karman  (1921)  boundary-layer 
integral,  except  that  xw,  here,  is  the  stress  as  it  appears  in  the  Cauchy  equation  (Slide 
9)  and  the  effect  of  the  body  force  (acceleration)  field,  fx,  is  explicitly  included. 

In  the  case  of  non-Newtonian  fluid  flow  and/or  in  the  presence  of  imposed  body 
(e.g.,  Lorentz)  forces,  it  is  important  to  include  (measure)  the  total  xw  not  just  the 
Newtonian  contribution. 

The  effect  of  injection/suction  velocity  at  the  wall,  vw,  is  noted.  Injection  can  be 
seen  to  decrease  shear  stress  at  the  wall.  However,  fluid  injected  with  zero  stream- 
wise  momentum  (as  assumed  here)  will  acquire  momentum  from  the  freestream, 
contributing  a  Reynolds  stress  and  increasing  drag. 

It  is  important  to  perform  integrated  drag-reduction  measurements.  In  the  case  of 
Lorentz,  or  other  spatially  inhomogeneous  force  fields,  shear-stress  at  the  wall  can 
be  driven  to  zero  locally,  even  as  there  may  be  an  overall  drag  increase.  Results 
from  some  measurements  we  reviewed  suffer  from  this  difficulty,  rendering  their 
assessment  problematical. 
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Momentum,  Body  Forces,  And  Stresses 


Momentum  transport  can  be  altered  by: 

*  Additional  body  forces 

-  Electromagnetic:  pf  -*  pf  +  je  x  B , 
where 

*  je  is  the  electric-current  density,  and 

*  B  is  the  magnetic  field 

*  Changing/adding-to  the  stress  field: 

Tcff  =  T(ViSC)  -  puV+  T(1)  +  T(2)  +  A 

where  changes  in  the  first  two  and/or  additional  stresses  can  be  the  result  of 
changes  in  the  fluid  rheological  properties,  e.g.,  via 

-  polymers 

-  surfactants 

-  microbubbles 


JASON  2000  Drag  Reduction 


Body  forces  and  additional  (non-Newtonian)  stress  contributions  can  oppose 
Reynolds-stress  contributions  and  lead  to  drag  reduction. 

In  the  case  of  polymers,  surfactants,  and  microbubbles  there  can  be  changes  to  the 
Newtonian  stress,  through  changes  in  the  shear  viscosity,  p,  In  comparing 
experiments,  it  is  important  to  distinguish  between  changes  in  drag  coefficients 
stemming  from  changes  in  (local)  Re  from  those  that  arise  from  changes  owing  to 
changes  in  rheological  behavior.  At  the  low  concentrations  that  high-MW  polymers 
are  typically  employed,  the  change  in  the  Newtonian  shear  viscosity  is  small.  This 
is  not  the  case  for  surfactants,  however,  that  are  used  at,  typically,  xlO  wppm  higher 
concentrations  than  high-MW  polymers. 

Care  should  be  exercised  in  assessing  experimental  results  to  make  sure  that  data 
are  correctly  parameterized  as  a  function  of  additive  concentration,  i.e.,  referenced 
to  the  correct  (modified)  Newtonian  shear  viscosity,  in  each  case. 
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Turbulent  Boundary  Layer  Drag 


Laminar-turbulent  transition  increases  wall  shear  stress  (drag) 

-  Turbulent/laminar  wall-shear  ratio  increases  with  increasing  Re 

-  Local  values  of  increased,  unsteady,  xw  can  be  partitioned  into 

•  a  higher  “background”  level,  with 

•  intermittent,  very  high  values  that 
represent  a  substantial  contribution 
to  the  mean  value. 

•  Latter  events  associated  with  lift  and 
bursting  of  “hairpin  vortices” 

-  ejection  of  low-speed  fluid  away 
from  wall  (large  Reynolds-stress 
contribution) 

-  Transient,  high-strain-rate  events. 

-  Frequency  ~U!  5  «  0.1-1  kHz  + 
at  U  ~  75  knots+ 


Image  by  Kamiadakis  (2000). 
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The  structure  of  near-wall  turbulent-boundary  layer  flows  and  identifiable  “events” 
responsible  for  momentum  transport  and  drag  are  subjects  of  continuing  research. 
“Hair-pin”  vortex-ejection  events,  which  bring  near-wall,  low-speed  fluid  away 
from  the  wall  have  the  correct  correlation  to  contribute  to  Reynolds  stress.  Such 
events  can  be  viewed  as  an  instability  of  the  counter-rotating,  streamwise  vortex 
structures,  which  are  spaced  by  V  =  X  /  v  «  100  and  centered  at/  =yux/v  =  20- 
30  (Blackwelder  and  Eckelmann  1979),  and  represent  intermittent,  energetic,  high- 
strain-rate  events. 

The  route  to  hair-pin  ejection  events  may  be  akin  to  that  of  the  Crow  instability 
(1970),  which  is  responsible  for  interrupting  the  lift-induced  counter-rotating 
trailing  vortex  pairs  (contrails)  by  pairing  interactions  behind  aircraft.  In  the  case  of 
the  wall-adjoining  counter-rotating  vortex  system,  the  self-induction  plus  wall- 
mirror  effects  will  result  in  a  lifting  of  vortex  lines,  as  indicated  in  the  computed 
image  by  Kamiadakis  in  the  frame. 

It  is  plausible  that  polymers  and  surfactants  act  by  resisting  the  transient  elongation 
(extensional  strain  rate)  associated  with  such  a  motion,  decreasing  the  frequency  of 
such  events  and  drag  thereby.  Such  resistance  to  extensional  strain  is  akin  to  the 
resistance  to  bending  of  magnetic-field  lines  in  MHD  flow. 
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Polymers  And  Surfactants 


Re^f 

Virk  (1975,/=  Cfl 


•  Flows  with  high-MW  additives, 

-  Polymers  (5-100  wppm) 

*  Polyethelene  oxide  (PEO,  e.g., 

C  Polyox  WSR301) 

‘43  *  Polyacrylamide 

gj  *  Guar  gum,  xantham  gum, . . . 

g  *  Carageenan  (seaweed  extract) 

-  Surfactants  (200  -  2000  wppm) 

*  Habon  G  (Ethoquad  T/l  3/50  with 

J*  sodium  salicylate  .,.) 

^2  *  Dobon  G 

t/i 

S  can  exhibit  high  drag  reduction 
(Toms  1948).  Demonstrated: 

-  [ACf/Cf]max~  60-80% 

•  Evidence  supports  improvement 
with  higher  Re' s,  with  limiting 
values  (as  long  as  polymers  remain 
intact) 

~  [ACf/Cf]max~  96% 
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Taking  the  ratio  of  the  Prandtl-Karman  friction-coefficient  relation  and  the 
maximum-drag-reduction  asymptote  suggested  by  Virk  (1975),  we  have, 

““(v/cTt,*  if  (vV)o  =>  max(Cf)pDl  «  0.04  (Cf)0  , 

for  large  Re's  (zero  subscripts  denote  “natural”  skin- friction  drag  coefficients). 
However,  at  fixed  flow  geometry,  the  maximum  Re  where  drag-reduction  benefits 
can  be  expected  will  be  limited  by  polymer  degradation  at  the  correspondingly  high 
strain  rates. 

A  somewhat  higher  drag-reduction  asymptote  has  been  reported  for  surfactants 
(Zakin  et  al.  1996). 
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Polymers:  Dynamics 


•  Polymers  are  chains  of  N*  104  -  105  monomer  ‘beads’  that  coil  to  produce  a 
cluster  of  size  (radius  of  gyration),  Rg  *  /  N°-6~  80  nm;  /  is  the  link  size. 

•  Phenomenological  (convected  Jeffreys,  or  Oldroyd-B  1950)  models  have 
introduced  convected  stress  fields  to  describe  viscoelastic  behavior 


-  Most  analyses  to  date  based  on  these  and  derivative  equations 
•  The  coiled  cluster  is  elastic  (wants  to  stay  quasi-spherical), 

-  with  a  free  energy  of  F  -  Vi  kh  R2 

•  kh=Ab  77 R2  is  the  effective  Hookean  (entropic)  spring  constant,  and 

•  Ab  the  Boltzman  constant. 

-  and  a  viscoelastic  relaxation  mechanism  (Zimm) 


Stokes  drag 


Ts 


dq 


q-Vu  =  - Knq  +  eT  =>  ®z  =  —  = 


kaT 


d  / 

elastic  force 


thermal 


r  \! 

excitation  \  7- 


Ys 


Zimm  frequency 


•  Zimm  (1956)  frequency,  coz,  provides  the  basic  relaxation  rate  that  must  be 
met  for  the  onset  of  polymer  viscoelastic  activity. 

-  coz*  105  rad/s  typical;  tz  ~  10'4  s  <=>  strain  rates  ~108  s_1  for  FastShip  BL 


JASON  2000  Drag  Reduction 


Polymer  molecular  weights  can  be  high  (  «  5-6x1 06  Daltons). 

The  convected- Jeffreys  and  Oldroyd-B  (1950)  models  are  discussed  in  Bird  et  al. 
(1987,  vol.  1).  Those  models  describe  coupled  convected  elastic  (tensor)  stress 
fields. 

Under  the  influence  of  thermal  excitations  (self-avoiding  random  walk),  long 
articulated  polymer  chains  coil  so  as  to  be  circumscribed  by  quasi-spherical 
boundaries.  The  coiled-chain  extent  is  measured  by  the  (rms)  radius  of  gyration, 

Rg.  Such  a  coiled  chain  will  deform,  when  subjected  to  hydrodynamic  (viscous) 
forces  in  a  strain-rate  field,  and  stretch  to  an  extent  R  >  Rg.  A  free  energy  of 
deformation,  F,  can  be  associated  with  such  elongations,  in  terms  of  an  effective 
(Hookean)  linear  (entropic)  spring  constant,  kh,  as  above. 

The  Zimm  frequency,  coz  =  KH/ys,  is  the  scaling  coefficient  of  the  ratio  of  the  left- 
hand  side  and  the  first  term  on  the  right-hand  side  of  the  equation  in  the  frame.  It  is 
a  function  of  the  temperature,  T,  the  fluid  viscosity,  p,  and  the  polymer-chain  size, 
Rg,  i.e.,  a  function  of  the  fluid  and  the  size  of  the  crumpled  polymer  chain,  but  not  a 
function  of  the  flow. 
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Polymers:  Drag  Reduction 


•  No  discernible  effect  in  laminar  (steady)  flow 

-  Only  small  changes  in  Newtonian  viscosity  at  concentrations  employed 

•  Modest,  if  any,  effect  on  laminar-turbulent  transition 

•  Drag-decrease  requires  wall  stress,  xw,  in  excess  of  critical  shear  stress, 
V 

-  Critical  wall  stress  independent  of  concentration 

•  For  x  >  x  : 

w  c 

-  Few  ppm  dramatically  reduce  drag 

*  Dilute  solution:  polymers  not  entangled 

-  Drag  reduction  increases  with  concentration  and  polymer  MW 

-  Experimental  evidence  of  maximum-drag-reduction  asymptote 

-  Profile  adjacent  to  viscous  sublayer  steepens 

-  u  V  decreases  in  steepcned-profilc  region 

*  u'  increases  slightly 

*  v'  decreases  slightly 

*  Reynolds-stress  decrease  is  not  the  result  of  decrease  in  turbulence  intensity 

JASON  2000  Drag  Reduction 


At  the  concentrations  employed  (a  few  wppm)  the  influence  of  polymers  on  the 
Newtonian  viscosity  of  the  fluid  is  small.  The  effective  viscosity  can  be  estimated 
from  the  Einstein  (1906)  equation  (cf  Landau  &  Lifshitz  6,  1959),  or  the  empirical 
Mark-Houwink  relation. 


lim 

Cp-»° 


/  \ 

Up  Vo) 


=  k,  m; 


In  which  and  c?g  are  empirical  constants  specific  to  the  polymer-solvent 
combination,  and  Mw  is  the  polymer  molecular  weight.  Virk  (1975)  tabulates  these 
coefficients  for  some  of  the  commonly  employed  polymers,  e.g.,  for  (polyethylene 
oxide)  PEO-water  and  Mw  in  units  of  106  Daltons,  Kyi  =  8.75*  10‘3  and  =  0.79. 

The  review  article  by  Virk  (1975)  can  be  consulted  for  experimental  evidence  of  the 
lack  of  discernible  influence  on  laminar-to-turbulent  transition  (Fig.  1),  a  critical 
wall  shear-stress,  xc,  that  is  independent  of  Re  (Fig.  2a),  and  polymer  concentration 
(Fig.  2b).  For  PEO,  Virk  reports, 


vnmy 


N/m2 


4.4  xlO6 


Documentation  for  the  drag-reduction  attributes  listed,  for  xw  >  xc,  can  also  be 
found  in  the  Virk  (1975)  review. 


19 


Polymers:  Drag-Reduction  Models  —  I 


Common  element:  Activation/extension  rate  o  Zimm  frequency 

W(^d)  U  n  1/2  u(X)  r  .  .  s 

- —  *  —  Re  >  -  >  coz  ,  for  X.  <  X  <  5 


dissipation  strain  rate  ^ 

Viscoelastic  layer - 


Wall 


inertial-eddy  strain  rate 


9  Y 

“X - 2T 

y.  x 


•  In  a  boundary  layer,  this  would  be  confined  to  a  viscoelastic  layer  between, 
yd  «  —  :  viscous/dissipation  distance 

Ux 

yz  «  —  :  distance  from  wall  where  Zimm  frequency  is  exceeded 
coz 
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Experimental  evidence  supports  the  notion  that  the  critical  shear  stress  is  related  to 
the  characteristic  time  (Zimm  frequency)  for  polymer  activation,  as  conjectured  by 
Lumley  (1973).  In  the  wall  region,  the  characteristic  length  and  time  scales  are, 

=  v  /  wT  and  f  /  uv  respectively  ( cf  Slide  14). 

Drag-reduction  onset  occurs  at  Rg  uj  v  «  8x1  O'3,  while  the  ratio  of  time  scales  is  tz 
ux2  /  v  «  2,  i.e.,  the  macromolecules  are  much  smaller  than  the  characteristic 
boundary-layer  (viscous)  length  scales,  but  have  a  (Zimm)  relaxation  time  (cf  Slide 
19)  that  is  comparable  to  the  viscous  time  scale  (Virk  1975). 

The  inequality  range  outlined  in  the  frame  derives  from  the  assumption  that 
Kolmogorov  scaling  can  be  used  in  the  vicinity  of  the  wall.  In  that  expression,  X,d, 
the  dissipation  length  scale,  is  assumed  approximately  equal  to  =  (v3/s)l/4,  the 
Kolmogorov  length  scale,  where,  here,  e  is  the  kinetic-energy  dissipation  (per  unit 
mass). 

That  polymers  are  not  expected  to  act  within  the  viscous  sublayer,  i.e.,  for/1  <  5,  or 
so,  was  also  postulated  by  Lumley  (1969,  1973).  The  resulting  viscoelastic  layer, 
foryd  <y  <yz,  then  defines  a  region  within  which  the  flow  is  not  damped  by 
viscosity  (Re  >  1)  and  where  strain  rates  are  high  enough  to  activate  the  polymers, 
defining  the  extent  where  local  momentum  transport  is  modified  (suppressed). 


20 


Polymers:  Drag-Reduction  Models  —  II 


•  Lumley  (1969, 1973): 

-  Identified  polymer-stretch  characteristic  time  as  criterion  for  activation 

-  Polymers  undergo  coil-stretch  transition  in  extensional  regions  of  the  flow 

-  Strongly  enhanced  viscosity  in  buffer  layer 
?  Turbulence  damping  <^>  drag  reduction 

*  Not  supported  by  subsequent  experiments  that  indicate  a  reduction  in  u-v 
velocity  correlation  magnitude,  even  as  u  and  v  fluctuation  levels  vary  only 
slightly 

*  de  Gennes  (1986)  and  Tabor  &  de  Gennes  (1986): 

-  Polymem  ‘elasticize’  turbulence 

-  Cascade  truncation  of  turbulence  cascade  when  F(h)  IV  ~  lA  p  u2(h  ) 

*  h :  elastization  length 

?  Drag-reduction  link  unclear 

*  Limited  utility  in  calculation 

JASON  2000  Drag  Reduction 


Lumley  correctly  theorized  that  polymers  act  beyond  the  viscous  sublayer  and  that 
drag-reduction  onset  is  associated  with  a  match  between  viscous  and  7.imm 
relaxation  time  scales  (Slide  21).  He  further  conjectured  that  the  net  effect  was 
equivalent  to  an  increase  in  effective  viscosity  and  a  kinetic  energy  damping,  with 
an  attendant  reduction  in  velocity  fluctuations  and  near- wall  Reynolds  stress.  The 
latter  conjecture  is  not  borne  out  by  experiment  (e.g.,  Virk  1975). 

Related  models  are  reviewed  by  Toonder  and  Nieuwstadt  (1999),  who  also  discuss 
an  idea  by  Joseph  et  ah  (1986)  and  Joseph  (1990),  who  argued  for  viscoelastic 
behavior  and  “shear  waves”,  and  conducted  some  experiments  to  look  for  them. 
Toonder  and  Nieuwstadt  conclude  that  this  idea  “. . .  is  only  vaguely  formulated  and 
has  not  been  experimentally  tested  or  theoretically  elaborated,” 

A  different  conceptual  framework  was  put  forth  by  de  Gennes  (1986)  and  Tabor  & 
de  Gennes  (1986),  who  based  their  proposal  of  polymer-flow  interaction  on  polymer 
deformation.  To  calculate  the  degree  of  extension,  they  invoke  Kolmogorov  scaling 
to  associate  a  unique  strain  rate  to  eddy  size.  They  hypothesized  that  the  turbulent 
cascade  is  “truncated”  at  scales  smaller  that  than  L  required  for  polymer  activation 
(see  frame),  at  equipartition  of  eddy  potential  (elastic)  and  kinetic  energy. 
Interestingly  (though  not  discussed  by  de  Gennes),  such  an  equipartition  implies  a 
competition  between  Reynolds  and  elastic  stresses,  with  a  concomitant  quenching 
of  turbulence  production  and  momentum  transport. 
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Polymers:  Elastodynamic  Equations  —  I 


Polymers:  Elastodynamic  Equations  —  II 


*  Elastization  <=>  drag  reduction  (Diamond,  Medvedev,  Nelson,  Gruzinov:  in 
preparation) 

-  Equations  do  not  predict  discernible  effect  for  laminar  pipe  flow 

-  Equations  predict  elastodynamic  waves  (akin  to  Alfven  waves) 

-  Elastization 

*  conversion  of  turbulence  kinetic  energy  to  elastic-wave  energy 

*  elastic/kinetic  energy  equipartition  (as  in  MHD) 

->  Elastic-tensor  stress  can  oppose  Reynolds  stress 

->rcff*-pwV  +  cpKnOX} 

*  Decreased  momentum  transport  for  y  <  yz 

-  Pro’s: 

*  Framework  for  implementing  de  Germe’s  intuition  on  elastization 

*  Recovers  momentum-flux  reduction  w/o  enhanced  damping 

*  Reveals  dependence  on  polymer-concentration,  cp 

-  Con’s: 

*  Restricted  to  Hookean  (linear)  springs 

*  Has  not  confronted  boundary-layer  structure/dynamics  (horseshoe  vortices,  etc.) 


JASON 2000  Drag  Redaction 
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Polymers:  Pro’s 

*  Small  amounts  needed  (10-50  ppm)  for  large  effects 

*  Proven  drag-reduction  capability 

-  Polymers  increased  peak  trans- Alaska 
pipeline  capacity  by  32,000  m3/day, 
at  10  ppm  levels 

•  injected  downstream  of  pumping  stations 
to  avoid  degradation 

•  No  loss  in  effectiveness  after  100’s  of  km 

•  Reduced  heat  transfer  (maintain  oil  viscos 

-  Used  in  storm  sewers  to  increase  peak  capacity 

*  Unaffected  by  sea  state  in  full-scale  (20-36%)  drag-reduction  trials 
(Highburton  1968). 

*  47%  drag  reduction  measured  at  64  knots  (Tuna  1971-1 974). 

*  Higher  drag  reduction  for  rough  walls 

*  Compatible  with  sea  water 
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Polymer  drag-reduction  figures  on  the  trans-Alaska  pipeline  system  and  on  storm 
sewers  are  from  Hoyt  (1990). 

Drag-reduction  data  cited  for  the  Highburton  (1968)  and  Tuna  (1971-74)  tests  are 
from  the  Ryan  (2000)  briefing.  The  latter  was  a  buoyant,  torpedo-like  vessel,  with 
25  m2  surface  area  and  is  significant  in  that  it  approaches  the  Reynolds  number  and 
strain-rate  regime  of  interest  in  the  Fast-Ship  context. 

That  higher  drag  reduction  is  observed  for  rough  walls  was  noted  by  Debrule  and 
Sabersky  (1974)  in  their  pipe-flow  experiments.  This  appears  to  be  corroborated  by 
experiments  in  progress  at  this  time  at  PSU/ARL  (Ryan  2000,  pvte.  comm.). 
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Polymers:  Con’s 

*  Some  logistical  complexity 

*  Unknown,  as  yet,  drag-reduction  etiology 

-  Extrapolation  to  bight’s  uncertain 

-  Little  guidance  for  optimization 

*  Rapid  mixing 

-  Higher  required  injection  rates 

-  Inhomogeneous  influence  along  streamwise  extent 

*  Polymers  break  (degrade)  at  high  strain  rates  (speed) 

—  When  local  strain  rate  exceeds  polymer  tensile  strength 

-  Stronger  polymers  can  be  made 

*  likely  to  degrade  more  slowly 

*  Injection  noise 


JASON  2000  Drag  Reduction 
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Surfactants  —  I 


•  Drag-reduction  benefits  recognized  by  Mysels  ei  al.  ( 1 949  patent,  +) 

•  Maximum  drag-reduction  asymptote  higher  (better)  than  for  polymers  (Zakin 
et  al.  1996) 

•  Dilute  solutions  (not  entangled)  at  concentrations  where  effective  drag- 
reducers 

-  «  x  l  o  wppm  high-M  W-polymer  concentrations 
>  Decrease  surface  tension:  relevant  to  microbubble-size  control 

•  Complex  molecules  with  hydro-phi lic/-phobic  ends  self-assemble  (reversibly) 
to  form: 


-  micelles, 

-  layered  structures, 

Depending  on  surfactant  geometry. 

-  vesicles,  and 

environment  salinity,  etc. 

-  worm-like  chains  v 

•  If  surfactant  assembly  is  broken  by  turbulence  (high  strain-rate  region), 

structure  self-heals,  in  regions  of  lower  strain-rate 
-  Used  in  closed-loop  heating  systems  to  reduce  pumping  costs 
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Surfactants  —  II 


*  Worm-like  chains: 

-  have  a  probability  distribution  of  contour  (chain)  lengths*  L 

where,  Pr(i)  *  e~L'k  £»  =  d,  cp'"'2  eE‘!{lk‘T> 

•  ~  4  nm  is  the  surfactant-worm  diameter 

•  <J>  ~  KH  surfactant-chain  volume  fraction 

•  £c~  0.37-0.40  eV  (  - 1 5-20  kBTt  at  room-71 )  scission  (tearing)  energy 

-  Persistence  length:  /p  »  16  nm 
—  Radius  of  gyration: 

Re  ~  *  (8nm)<p1/4  eE<m*T) 

-  Number-density  of  surfactant  polymers: 

cp*(2.5xl020cm-3)<p1/2e-£«™ 

-  Non-Newtonian  (viscoelastic)  stress  contribution  (as  for  high-MW  polymers): 

•  t<Po1)  =  cp  kh  Q,  with  kh  =  kBT  i  i?g2  the  Hookean  spring  constant 
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Polymer-Surfactant  Comparison 


PEO  (polymer) 

C14TABr  (surfactant) 

Drag  reduction 

~  4-fold 

-  3 -fold 

Persistence  length:  /p 

-  1.5  nm 

~  16  mn 

Contour  (chain)  length:  L 

-  104nm 

~  350  nm 

Segments:  N=  L  / 1? 

-6500 

-20-25 

Radius  of  gyration:  i?g 

~  300  nm 

~  75  nm 

Polymer  concentration:  Cp 

~  1014  cm'3 

-  3xl014  cm*3 
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Microbubbles 
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Microbubble  injection  results  in 
comparable  drag  reductions 
(McCormick  and  Bhattacharya 
1973) 

-  [ACf  /  Cf]nmx  =  60  -  80% 

-  Similarities  in  phenomenology  to 
polymers/surfactants . 

-  However: 

*  Comparable  gas-volume  fractions 
(-80%)  found  to  be  required  for 
reduction  at  such  levels 

*  A  consequent  reduction  in  the 
local  two-phase  fluid  density,  p, 
implies  a  comparable  reduction  in 
local  Reynolds  stress:  -  pu  V 

High  pumping  costs 

-  Comparable  (on  paper)  with 
anticipated  drag-reduction  benefits 
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Microbubble  Dynamics  —  I 


•  Microbubble  size  distribution  adjusts  to  local  flow: 

-  Breakup:  Balance  between  capilarity  and  shear-stress/pressure  forces  (Hinze  1955) 

-  Coalescence:  When  film  between  adjacent  bubles  drains  to  critical  thickness,  /?, 
and  ruptures  (Thomas  1981) 

0.7  (a/p)3/5  s"2/5  ~  xw3/5  <  <  2.4  (o2/7/pp)3/5  zA!A-  xw*3/8 

with  a  the  surface  tension  and  s  -  wT3  /  62  the  BL  kinetic-energy  dissipation  rate 

•  Microbubbles,  at  void-fractions  of  interest  (comparable  to  %  drag  reduction), 
substantially  lower  sound  speed 

-  Compressibility /turbulence-damping  effects 

-  Conversion  of  TKE  to  acoustic  radiation;  significant  modification  of  BL  acoustics 

-  “Supersonic”  boundary  layer 
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Microbubble  Dynamics  —  II 


•  Microbubbles  can  impart  elastization  to  the  flow  with  phenomena  that  parallel 
polymer/surfactant  rheological  consequences  and  dynamics 

-  Smaller  bubbles  more  effective 

*  p-reduction  from  void  fraction  not  the  whole  story 

*  Reduction  in  required  volume  fraction 

-  Wall  spectral  data  from  pressure  and  hot  films  show  loss  of  high/fluctuations 

>  f  P  s4)3/5  /  v3/4  >  1  Harder  as  U  f 

-  Polymers  and  microbubbles  both  reduce  turbulence  statistics  skewness 

*  Decrease  of  high-amplitude  events 

*  Electrolytic  microbubble  generation? 

-  First  by  McCormick  and  Bhattacharya  (1973) 

-  Generate  bubbles  at/near  wall,  where  needed 

*  Microbubbles  generated  within  y+  <  10. 

-  Gas  volume  flow  rate  estimated  to  be  much  lower  (1 0  4 -1 0  s)  than  for  porous-wall  injection 

-  Pulsed/adaptive? 

-  Attractive,  but  net  power  gain  may  be  difficult: 

*  Low  seawater  electrical  conductivity,  ae. 
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Air  Films 


•  Air  films  isolate  high-Vu  near-wall  region  from  water. 

-  Reduced  tw  =  [|i(d«%)]w  ,  since  pair  *  |iwatcr  /  80 

•  Can  be  more  effective  than  polymers/microbubbles: 

-  If  stable 

•  at  high  Re' s 

•  Rayleigh-Taylor  instability 

•  Interfacial  waves 

•  Freestream  turbulence 

-  If  pumping  costs  can  be  recovered 

-  JSR-99-215  proposal 

•  Exploit  high-pressure  ratio  at  TE 
to  drive  supersonic  flow  out  of 
convergent-divergent  nozzle 

•  Graceful  “failure” 

-  Air  film  converted  to  microbubble  layer 
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E&M  Turbulence  Control  (EMTC)  —  I 


*  Exploit  conductivity  of  sea  water 

*  Lorenz  force: 

-  Permanent/electro-magnets  in  wall 

-  f  =  jexB*0eExB  may  be  streamwise,  spanwise,  or  wall-normal 

-  Wall  electrodes  can  drive  currents  L  B 

*  Does  it  work? 

-  Sometimes  drag  increase,  sometimes  drag  decrease 

*  Cannot  predict  which 

*  Empirical  —  no  theory  for  turbulent  flow 

*  Depends  on  (many)  specific  parameters 

-  Both  stationary  and  travelling- wave  E&M  fields 

*  Experiments  (Brown):  up  to  ~1 5%  reduction  reported 

*  Simulations  (Kim  &  Kamiadakis):  up  to  -35%  reported  (low-/moderate-Re) 

-  Unknown  scaling  with  Re 
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E&M  Turbulence  Control  fEMTQ  —  II 


•  Efficiency:  r)EM  =  ~  °e^  ^ 

*EM  P», 

-  |B|  limited  to  ~  0.5  T  (5  kgauss) 

-  Low  conductivity  of  seawater  (oc  -  4.3  mhos/m  —  PC:  Siemens/m) 

=>  Large  E  =^>  large  resistive  loss 

-  CFD  (Kim,  Kamiadakis)  indicates  r|EM  »  10  3. 

•  Can  it  be  improved  to  be  useful? 

-  Closed-loop  control:  apply  force  only  where/when  useful 

-  Adaptive/leaming  (neural-net)  algorithm? 

-  CFD  indicates  improved  efficiency  to  rjEM  ~  1/6  (Kim) 

-  Unknown  limits  and  whether  r|EM  >  1  is  possible. 

•  Practical  issues 

-  Skin  covered  with  electromagnets/electrodcs 

*  1 08  —  1 010  electrodes  for  full-scale  vessel 

-  Low  voltage  =>  high-current  supply 
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Active-Feedback  Control 


*  Pro's: 

-  Effectiveness  demonstrated  (simulation)  in  low-j?e  turbulent  boundary  flows 

-  Potential  with  emerging  MEMs  technology 

*  Con's: 

-  Actuation  and  efficacy  unknown  at  high  Re*  s 

-  Effective,  realizable,  survivable  actuation  as  yet  undefined 

-  Unknown  efficiency  (Power  saved  vs.  expended) 
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Conclusions 


*  Required  FastShip  *3-4  drag-reduction  levels  appear  to  be  within  reach. 

*  However,  real  requirement  {net  power  reduction,  for  given  U)  is  harder: 


Reduced  CD  ■ 


CD  + 

X' 


pU  AH 


<  /mC.  with  /DR  <  0.25-0.30 


“Natural”  CD 
■  Normalized  drag-reduction  power-cost 


•  PDR  is  the  power,  from  the  primary  power  plant,  to  effect  drag-reduction  (in  the 
case  of  EMTC,  possibly  with  its  own  efficiency  penalties) 

*  Care  required  in  designing/interpreting  experiments/simulations  to  account  for  proper 
(total)  shear  stress  at  the  wall  (drag)  and  PDRi  drag-reduction  power-expended 

•  Generalized  von  Karmdn  momentum  integral 

•  g  ( du/dy ^  may  not  include  necessary  contributions  from  viscoelastic  stresses  for 
polymers/surfactants,  or  Lorenz  forces  (in  case  of  EMTC) 

*  Proposed  unified  theoretical  framework  to  account  for  behavior  of 

•  high-MW  polymers  and 

•  Surfactants 

•  High  community  interest:  35,108  “drag-reduction”  hits  on  Netscape 
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Recommendations  —  I 


*  Solve  turbulence  problem 

-  Was  not  accomplished  this  My 

•  Polymers: 

-  Proposed  elastodynamic  equations 

*  Boundary  conditions?  Steady-flow  response? 

*  Linear-stability  analysis  to  verify/validate  onset  phenomena 

*  Computer  simulations  (DNS)  to  investigate  implied  drag-reduction  mechanisms 

*  Validation  experiments  (e.g.,  pipe  flow)  of  predicted  Alfvdn-like  waves 

'  -  Dispersion  relation;  co(k)  -  -  2icoz  ±  [<ep  k#  Qq  f  p)  k2  -  4&z2]m 

-  Influence  on  turbulence;  k  >  =  2&z  /  (cp  ku  Q&  k1 1  p)m  attenuated. 

*  Experimental  investigation  of  modifications  to  bursting  events  and  near-wall  flow 
structure 

-  Extend  elastodynamic  framework 

*  Lift  limitation  of  linear  response  (infinite  extensibility) 

-  Linear  response  does  not  capture  measured  non-Newtonian  (shear-rate-dependent) 
shear  viscosity 

-  Undertake  hybrid  (continuum/partiele)  simulation  of  non-linear  elastodynamic 
“dumbell- fluid” 

-  Contribution  of  polymer  stiffening  at  large  extensions  (prior  to  breaking) 

-  Likely  to  be  important  in  impeding  Mgh-strain-rate  events  (hair-pin  vortices) 

-  “Designer  polymers”;  Prescribed  stress-strain  relation? 

JASON 2000  Drag  Reduction 


37 


Recommendations  —  II 


•  Surfactants 

-  Exploit  self-assembly/-healing  to  tailor  non-Newtonian  behavior  where  needed 

•  Flow-adaptive  response 

-  Use  in  closed-loop  heat-exchange  systems  to  reduce  power/noise 

•  Microbubbles 

-  Explore  flow-elastization  idea 

-  Develop  unified  theoretical  framework  for 

•  polymers 

•  surfactants 

•  microbubbles 

•  Vapor  cavities  and  air  films 

-  Stability  issues  at  high  Re's 

—  Extend  present  computational  capabilities  to  accommodate  dynamic-pressure 
variations  in  cavitating/ventilated  regions 

-  Study  air-film  injection-power-recovery  technologies 

-  Air-film  — »  microbubble  transition 
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Recommendations  —  III 


*  If  hydrofoil  based,  Fast-Ship  DR  technology  likely  to  be  composite: 

-  Optimized  to  different  flow/pressure  regimes 


-  Revisit  system  optimization,  applying  potential  DR  benefits  to  various  parts  to 
explore  changes  in  shape,  structure,  and  consequent  system  L  ID . 
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